The p meson decay constant using a tadpole-improved action 
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The p meson decay constant and the associated renormalization factor are computed in the quenched approxi- 
mation on coarse lattices using a tadpole- improved action which is corrected at the classical level to 0{a^). The 
improvement is displayed by comparing to Wilson action calculations. 
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1. INTRODUCTION 

By improving a lattice action, one hopes that 
calculations at a fixed lattice spacing will be 
closer to their continuum values than with the 
original action, and similarly that calculations of 
a fixed accuracy can be performed on coarser lat- 
tices than those required, for example, by the 
Wilson action. For the present work, we con- 
sider an action which has classical errors begin- 
ning only at 0{a?), and which is also tadpole- 
improved This action has been used to com- 
pute the spectrum of light hadron masses |^,||, 
and it was found that hadron mass ratios arc con- 
siderably improved compared to Wilson action re- 
sults at comparable lattice spacings. In this work 
we address the question of whether matrix ele- 
ments are similarly improved by considering the 
vector current renormalization factor and vector 
meson decay constant. 



2. IMPROVED ACTION 

We use an improved gauge field action which 
involves a sum over 1x2 plaquettes as well as 
1x1 plaquettes and an improved fermion ac- 
tion which contains next-nearest-neighbour inter- 
actions llsll, 
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The tadpole factorg is Uq = {^ReTiUpiY 

The presence of three timesteps in a single term 
of the action introduces artifacts]^, such as oscil- 
lations in meson mass functions near the source. 
However on our lattices, the plateau itself is not 
altered. 

3. METHOD 

In the continuum, the p meson decay constant 
fp is defined as follows: 

(0|Fp|p)cont = fp^mlcp . 



(2) 



On a lattice, matrix elements of the vector current 
get renormalized by a factor Zy, 



{f\Vp\t),ont=Zv{9^){f\V^\t) 



(3) 



where the dots represent finite lattice spacing ef- 
fects. For the improved action, the 0{a, a^) terms 
vanish identically and terms proportional to the 
coupling are kept small by the tadpole factor. By 
considering the ratio of the local vector current 
Vi^{x) = 4'{x)"ipip{x) to the conserved vector cur- 
rent {x)^ we can compute the renormalization 
factor, 



Zv{9') = {fK\^l{fK\t) + 



(4) 



We will consider a variety of choices for the initial 
and final states in this ratio, and the differences 
in the resulting values for Zy will reflect the finite 
lattice spacing effects. 



4. CALCULATION 

We have performed calculations at two values 
of (3 with the improved action which correspond 
to lattice spacings of 0.4fm and 0.27fm as de- 
rived from the string tension. For comparison, 
we have also used the Wilson action to calculate 
at the same lattice spacings. For the gauge fields 
we used pseudo-heatbath updating with periodic 
boundary conditions; the first 4000 sweeps were 
discarded and then 250 sweeps (200 sweeps) were 
discarded between each pair of improved (Wilson) 
configurations that was kept. 

A stabilized biconjugate gradient algorithm 
was used for the fermion matrix inversion at three 
choices of the hopping parameter. Due to the 
modest number of lattice sites, Dirichlet bound- 
ary conditions were used for fermions propagat- 
ing in the time direction but periodic boundaries 
were used in all spacial directions. The source was 
placed two timesteps away from the boundary. 

Table 1 

Some details of the simulations. Njj is the num- 
ber of gauge field configurations and a^t is the 
lattice spacing derived from the string tension. 



Lattice 


Nu 


13 ast [fm] K. 






Improved Action 


6^x12 


100 


6.25 0.4 0.162, 0.168, 0.174 


8^x14 


50 


6.8 0.27 0.150, 0.154, 0.158 






Wilson Action 


6^x12 


75 


4.5 0.4 0.189, 0.201, 0.213 


8^x14 


30 


5.5 0.27 0.164, 0.172, 0.180 



5. RESULTS 

Figure |l| shows our determination of Zy from 
four separate choices of the states in Eq. (|4|), 

< tt\V\'k >, < p\V\p >, < N\V\N >, < Q\V\p > . 

The extrapolation to the chiral limit was done lin- 
early in toJ. The results of Wilson calculations 
by other groups are also shown. Notice that for 



N* 0.50 



- 


i 


1 1 1 1 1 1 1 1 

l-«H 


1 1 1 1 1 1 

- 

* ■ 




o 






■ A Ref. 


[7] 






■ Ref. 


[8] 






; □ Ref. 


[9] 







0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

M Q 
p 

Figure 1. The renormalization constant in the 
chiral limit. Solid (open) symbols are from the 
improved (Wilson) action. The lowest lying point 
at each Mpa is from the decay amplitude. 



both actions, Zy at fixed MpU is approximately 
independent of the amplitude used, except for the 
decay amplitude which predicts a smaller value 
for Zy- For fixed lattice spacing, the improved 
action does give improved values of Zy', in fact 
the improved action matches the Wilson predic- 
tions for Zy when the lattice spacing is larger by 
a factor of about three. Figure || shows that the 
dependence on Af^ is also quite similar for the 
two actions. 
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Figure 2. The renormalization constant as a func- 
tion of . Solid symbols are from the improved 
action at /? = 6.8; open stars and diamonds are 
Wilson at (3 = 6.0; open squares are Wilson at 
P = 5.85. 
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Figure 3. The lattice decay constant. The dashed 
Une is the experimental value; solid symbols are 
from the improved action; open symbols are Wil- 
son calculations. All Wilson quarks are normal- 
ized by \/2k. 
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Figure 4. The renormalized decay constant. The 
dashed line is the experimental value; solid sym- 
bols are from the improved action; open triangles 
and open squares are from the Wilson action. All 
other symbols are Wilson calculations with Zpat ■ 



Figure |3| shows the lattice decay constant (i.e. 
without Zv)- Again, we find that the improved 
action does give improved results when compared 
to Wilson calculations at the same lattice spacing, 
and that the improved prediction is similar to the 
Wilson prediction at 1 /3 of the lattice spacing. 

Our best determination of the decay constant 
comes from the product of Zy and l/(Zy/p), and 
this is shown in Figure ^. The cluster of data 
points which lie above the experimental result 
rely on tadpole-improved perturbation theory to 
estimate Zy- Notice that the completely non- 
perturbative determinations of l/fp are quite flat 
over a substantial range of lattice spacing, and 
that for both the Wilson and improved actions 
the result is near experiment. 

In summary, we have found that the C(a^) 
classical and tadpole improvements of our action 
manifest themselves in vector current matrix ele- 
ments in essentially the same way as in masses . 
Improved action results for l/fp and Zy are con- 
sistent with Wilson action computations done at 
a lattice spacing about 1/3 the size. 
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